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Abstract 

In a variety of models of neutrino masses and mixings the lighter top squark decays into 
competing R - parity violating and R - parity conserving channels. Using Pythia we 
have estimated in a model independent way the minimum value of P = BR(ti ^ ex?) x 
BR(ti ifb), where li = and fi^, corresponding to an observable signal involving 
the final state 11 + jets +fjT (carried by the neutrinos from the x? decay) at Tevatron 
Run II. For the kinematical cuts designed in this paper P depends on m~ only. We 
then compute P for representative choices of the model parameters constrained by the 
oscillation data and find that over a significant region of the allowed parameter space 
P is indeed larger than Pmin- This signal is complementary to the dilepton + dijet 
signal studied in several earlier experimental and phenomenological analyses and may 
be observed even if BR(ti ifh) is an order of magnitude smaller than BR(ti cxi). 
The invariant mass distribution of the hardest lepton and the hardest jet may determine 
mr and reveal the lepton number violating nature of the underlying interaction. The 
invariant mass distribution of the two lowest energy jets may determine mrro. 
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1 Introduction 



Neutrino oscillations in different experiments have [Ij established that the neutrinos have 
tiny masses, several orders of magnitude smaller than any other fermion mass in the Standard 
Model ( SM ) with massless neutrinos. Massive Dirac neutrinos can be accommodated in 
the SM if right handed neutrinos are introduced as SU(2) singlets. But the corresponding 
Yukawa couplings must be unnaturally small. There are several more aesthetic mechanisms 
of introducing neutrino masses. We shall briefly review below two popular approaches based 
on supersymmetry (SUSY) [2] and their possible impacts on the high priority program of 
SUSY search at high energy accelerators - on the current and future experiments at Tevatron 
Run II in particular. 

The see-saw mechanism [3] in a grand unified theory(GUT) with or without super- 
symmetry (SUSY) [2j, offers a natural explanation of small neutrino masses provided the 
neutrinos are Majorana fermions. One need not fine tune the Dirac masses to unnaturally 
small magnitudes. Instead a typical neutrino Dirac (Majorana) mass in this model is as- 
sumed to be of the order of the electroweak scale(GUT scale (M^)). The physical neutrino 
masses turn out to be proportional to the ratio of these scales of widely different magnitudes 
and are, therefore, naturally suppressed. 

The observation of neutrinoless double beta decay [5] will provide a strong indirect 
evidence in favour of the Majorana neutrinos. Another hall mark of any GUT is the proton 
decay [1] which has not been observed so far. However, all non-supersymmetric GUTs suffer 
from the naturalness problem [2] which destabilizes the mass of the higgs boson essentially 
due to the same large hierarchy of the two mass scales responsible for the see-saw. 

A super symmetric GUT (SUSYGUT) cures the hierarchy problem provided the masses 
of the sparticles (the supersymmetric partners of the SM particles) are O (1 TeV). Thus 
the exciting program of sparticle searches and the reconstruction of their masses at the on 
going ( Tevatron Run II) and the upcoming (the Large Hadron Collider (LHC) or the Inter- 
national Linear Collider (ILC)) accelerator experiments have the potential of testing SUSY. 
Furthermore, in simple grand desert type models SUSY indeed facilitates the unification of 
the three couplings of the SM at a scale compatible with the current constraints from proton 
decay. 

It should, however, be emphasized that in the most general framework with a chosen 
GUT group the neutrino masses and mixing angles involve many unknown free parameters 
(e.g, the elements of the Dirac and Majorana mass matrices). Collider experiments provides 
very little information on this sector. On the other hand Neutrino data alone cannot fully 
determine these parameters unless additional assumptions are introduced to simplify the 
neutrino mass matrix [B]. Such assumptions involve the minimal choice of higgs multiplet 
in a GUT, imposition of additional discrete symmetries etc. Thus the observation of both 
neutrinoless double beta decay and proton decay along with the discovery of sparticles with 
masses at the TeV scale may at best be regarded as a circumstantial evidence of an underlying 
SUSYGUT. There is no simple way of relating the measured sparticle masses with the physics 
of the neutrino sector. 
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Although theoretically the idea of unification of the couplings is rather appealing there is 
no compelling experimental evidence in favour of it. More importantly if neutrinoless double 
beta decay is observed but the proton decay remain illusive, the case for an alternative 
theory would be strengthened. Such a theory is provided by the R-parity violating (RPV) 
super symmetry, where R-parity is a discrete symmetry under which the particles and the 
sparticles transform differently [7j. It should be noted that the minimal supersymmetric 
extension of the SM (MSSM) naturally contains a R - parity conserving (RPC) as well as a 
RPV sector [2j . But the couplings in the latter sector violate both lepton number and baryon 
number resulting in catastrophic proton decays. Thus one option is to impose R-parity as a 
symmetry and eliminate all RPV couplings. This model is generally referred to as the RPC 
MSSM. 

But there is an alternative. If either baryon number or lepton number conservation 
but not both is required by imposing appropriate discrete symmetries then the catastrophic 
proton decay is suppressed. Such models are usually referred to as the RPV MSSM. 

The lepton number violating version of the RPV MSSM is more appealing since it natu- 
rally leads to Majorana masses of the neutrinos [3 El E] and neutrinoless double beta decay 
0.1 

More importantly, the observables in the neutrino sector in this depend not only on the 
RPV parameters but also on the RPC ones (including the sparticle masses). Thus the precise 
determination of the neutrino masses and mixing angles in neutrino oscillations and related 
experiments on the one hand and the measurement of sparticle masses and branching ratios 
(BRs) at accelerator experiments on the other, can indeed test this model quantitatively. In 
addition the collider signatures of this model are quite distinct from that of the RPC model 
(see below). In this paper our focus will be on a novel signature of a RPV model of v mass 
which can be probed at Run II of the Tevatron collider. 

In the RPC MSSM the lightest supersymmetric particle (LSP) is necessarily weakly 
interacting, stable and, as a consequence, a carrier of missing transverse energy [fir)- This 
^ is a hallmark of RPC SUSY. In contrast the LSP is necessarily unstable and decays into 
RPV channels, violating lepton number in the model under consideration. In addition other 
sparticles can also decay into lepton number violating modes providing very novel collider 
signatures. The BRs of the latter decays, however, depends on the relative magnitudes 
of their widths and that of the competing RPC channels. Thanks to the LSP decay the 
multiplicity of particles in any event is usually much larger compared to the corresponding 
event in RPC SUSY containing the stable LSP. Moreover, since the LSP is not a carrier of 
frr the reconstruction of sparticle masses appears to be less problematic in a RPV model. 
We shall consider a few examples of such reconstructions later. 

Apparently the stringent constraints from the neutrino data (discussed below) implies 
that the RPV couplings must be highly suppressed. This in turn suggests that the branching 
ratios (BRs) of the lepton number violating decays of sparticles other than the LSP will be 

■^The RPV MSSM can be accommodated in a GUT by introducing, e.g., non-renormalizable higher 
dimensional operators; see the concluding section for further discussions and references. 



2 



very small compared to the competing RPC decays. One notable exception, however, is the 
direct RPV decay of the lighter top squark (ti) [lOl [TTl [121 [13], if it happens to be the next 
lightest supersymmetric particle(NLSP) while the lightest neutralino {xi) is the LSP. The 
ti-NLSP assumption is theoretically well-motivated due to potentially large mixing effects 
in the top squark mass matrix [2]. In this case the RPC decays of ti occur only in higher 
orders of perturbation theory and are also suppressed (to be elaborated later). Thus they 
can naturally compete with the RPV decays even if the latter modes have highly suppressed 
widths as indicated by the neutrino data. Thus the competition among different decay modes 
of the lighter top squark, which may be the only strongly interacting sparticle within the 
striking range of Tevatron Run II, is a hallmark of RPV models of neutrino mass [T3] . 

No signal of either RPC or RPV has been observed so far. Thus constraints on the 
RPV parameters O [H] have been obtained from the experimental data. As in the case of a 
SUSYGUT, the most general RPV framework has too many parameters. Thus one usually 
considers some benchmark scenarios, each consisting of a minimal set of RPV parameters at 
the weak scale [151 [IS]! which can produce an acceptable neutrino mass matrix consistent 
with the oscillation data. 

Among the examples in ref. ^T^, we have focused on a specific model with three trilinear 
couplings A ( where i = 1,2,3 is the lepton generation index ) which can trigger ti decays 
and three bilinear RPV parameters Hi at the weak scale. With the above couplings the 
neutrino masses turn out to be proportional to m^. A brief review of this model along with 
notations used here may be found in section 2 of [T7] . 

The stringent upper bounds [ISl [IZ] on the trilinear (bilinear) couplings from neutrino 
data are ~ 10~^ ( ~ 10~'^GeV ). Thus almost all collider signatures arising from these 
couplings except the LSP decay are expected to be unobservable. As already mentioned, a 
notable exception could be the direct RPV decay of the ti - NLSP via a A^33 type coupling 
[121 [13] into a 6-quark and a charged lepton. 

tl - Itb, (1) 

where i = e, /i or r. 

This is so because the competing RPC decay modes in this case are i) the loop induced 
decay [T^ 

h cx? (2) 

and ii) the four body decay [12] 

tl - &X?/f (3) 

where f — f refers to a pair of light fermions. The last two decays occur only in higher 
orders of perturbation theory and, consequently, their widths are also highly suppressed. It 
may be recalled that for relatively large values of tan f5 the loop induced decay overwhelms 
both the RPV decay [13] and the four body decay [HI [201 [21] • 

In the earlier experimental and phenomenological analyses signals from the pair pro- 
duction of ti-tl [m 1121 US] followed by RPV decays of both were considered. The model 
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independent minimum observable branching ratio (MOBR) of the channel ti — > e^b for Teva- 
tron Run II was also estimated as a function of the lighter top squark mass {mj^) [13] by 
considering the dilepton + dijet signal. We shall review these limits and their consequences 
in a later section. Moreover, it was also demonstrated that by reconstructing the two lep- 
ton - jet invariant masses the lepton number violating nature of the decay can be directly 
established [13] for a variety of mj^ . 

In this work we shall concentrate on a new signal which arises if one member of the ti-tl 
pair decays via the RPV channel into an electron or a muon while the other decays via the 
loop induced channel (Eq.2). The second decay is assumed to be followed by the LSP decay 
via the modes 

x'^-^u,bb ; z = 1,2,3 (4) 

The modes in Eq.4 indeed occur with a combined BR of 100% in all models provided 
the LSP mass m~o is smaller than mw (the W boson mass) so that its decays into W, 
Z or t are kinematically forbidden. In addition it holds to a high degree of accuracy for 
any m~o if the higgsino components of the LSP are highly suppressed(e.g., if it is almost a 
pure bino) and it is lighter than nit B Thus the signal consists of a very energetic lepton 
accompanied by several jets (we do not employ flavor tagging) plus a moderate amount 
of missing energy carried by the neutrinos. We have generated the signal events by using 
Pythia(version 6206) [22] • The background events, discussed in detail in the next section, 
are generated either directly by Pythia or by interfacing Pythia and CalcHEP (version 2.3.7) 
[23j . We then introduce kinematical cuts which optimally suppress the backgrounds. Finally 
we define the model independent product branching ratio (PER) 

P = BR{ti l+b) X BR{tl cXi), (5) 

where BR(ti l^b) = BR{ti e^b) + BR{ti fi^b) and Xi is assumed to decay 
into a pair of jets and missing energy (Eq.4) with 100 % BR. Our simulations estimate the 
minimum value of P {Pmin) as a function of mj^ corresponding to an observable signal. These 
estimates are for an integrated luminosity of 9 fb~^ and S/ \fB = 5 where S and B are the 
total number of signal and background events. We have used the top squark pair production 
cross section as given by QCD [21] as an input. 

We have also examined the invariant mass of the most energetic lepton and jet in the 
signal for various mj^ . The resulting distribution peaks around the input value of mj^ . This 
illustrates that the combinatorial backgrounds are not very damaging. The peak, if observed, 
will unambiguously demonstrate the lepton number violating nature of the underlying inter- 
action. 

Our estimates also establish that the signals proposed in [13] and that proposed in this 
paper are complimentary. While former happens to be the most promising search channel if 
the RPV BR dominates over that of the loop induced decay, the latter can potentially reveal 

"^see section 3 for further details. 
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the presence of RPV interactions even if the loop decay overwhelms the RPV decays. These 
points will be further illustrated in the subsequent sections with numerical examples. 

We then turn our attention to some specific models. We compute P in these models for 
the entire RPV parameter space allowed by the oscillation data and compare the results with 
Pmin estimated by our Monte Carlo. We find that for top squark masses within the striking 
range of the Tevatron, a large region of the allowed parameter space (APS) yields P greater 
than the estimated Pmin- We have also checked that in significant regions of this parameter 
space the BR(ti — Ifh) is in fact smaller than the estimated MOBR in or the updated 
value given in section 3. Thus the signal introduced in this paper is indeed complementary 
to the ones studied earlier and may turn out to be the main discovery channel of ti. 

The plan of the paper is as follows. In section 2 we present our Monte Carlo studies 
of the signal and the backgrounds. This is followed by the main results of this paper: the 
model independent estimates of Pmin as a function of mj^ . A comparison of the viability of 
this signal vis a vis that of [13] is also included. In section 3 we compute the parameter P 
in some specific models and revise the estimated MOBR in [13j for representative parameter 
spaces consistent with neutrino oscillation data. This study establishes that the two signals 
are indeed complementary. Our conclusions and future outlooks are summarized in the last 
section. 



2 The Signal and the SM Background 

We have simulated the signal described in the introduction by generating 10^ events using 
Pythia. The most energetic I — h pair in the final state comes most of the time from the 
direct decay of ti. The kinematics of this pair is, therefore, independent of fn-o. Using this 
feature we have shown that it is possible to choose the optimal kinematical cuts in such a 
way that the efficiency of these cuts is practically independent of rn~o^. The size of our signal 
is completely determined by i) the cross section of ti pair production as given by QCD [21] 
and the efficiency which depend only on mj^ among the SUSY parameters and ii) the model 
independent input value of P (Eq.5). We have considered the following backgrounds: 

1. tt, 2.W+W-, 'i.WZ, 4.WH, 5.bb, Q.th + ib, 7.W + 2] 

Backgrounds 1 to 5 have been simulated by Pythia. Backgrounds 6 and 7 have been 
generated by CalcHEP at the parton level. Subsequently initial and final state radiation, 
hadronization, decay and jet formation have been implemented by interfacing with Pythia. 
All cross sections are calculated by CalcHEP in the leading order using the CTEQ6M par- 
ton density functions [25] using the four flavour scheme. The next to leading order(NLO) 
corrections would modify both the signal and the backgrounds by the appropriate K-factors. 
For example, the recent QCD prediction for the K-factor for the ti cross section, the most 
dominant background (see Table [1]), is 0.94 to 1.52 depending on the choice of the renor- 
malization scale in the leading order cross section |26j. The corresponding number for the 
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signal is ~ 1.3 [27]. Since our main result( the estimated Pmin) is based on the ratio S/^/B, 
we beheve that neglecting the NLO corrections would not change our conclusions drastically. 

For the background from Wbb, a has been computed with nominal cuts of Pt > 3 GeV 
and \rj\ < 4.0 on the parton jets. This is to eliminate the soft and collinear processes which 
are important for the NLO calculation. Here our main aim is to generate events with high 
Pt, central jets which can contribute to the background surviving the cuts listed below. 
We have checked that reasonable variations of these nominal cuts do not influence the final 
results. 

The backgrounds from Wss and Wcc have not been simulated separately since flavour 
tagging is not included in our selection criteria and for our cuts these contributions are 
expected to be similar to that of Wbb. Similarly the Wdd and Wuu backgrounds computed 
with the same nominal cuts as above are practically identical and only one of them has been 
simulated. Unless otherwise mentioned contributions of W's of both signs are included in 
Table [1]. 

We present the kinematical cuts and their efficiencies for the signal and the backgrounds 
in Tables [1] and [2] for mj^= 180 GeV. The cuts (Ci — C3) are as follows: 

1. Ci:Only events having an isolated lepton (e or /i of either charge) with AR{l,j) > 0.5 
, \rii\ < 2.5 and Pt > 105 GeV are accepted. Here AR{m.n) = ^/A(|)^~+~Arf , m,n stand 
for either a lepton(l) or a jet(j). 

2. C2:The number of jets in an accepted event is required to be rij > 2, where, jets are 
selected by the toy calorimeter of Pythia if Ej, > 12 GeV, \rij\ < 2.4 and AR{ji,j2) > 
0.5. 

3. CaiEvents with the invariant mass of the two highest Pt jets lying between 100 GeV 
< Mj^j^ < 70 GeV are rejected. 

The cuts are applied in the order shown in the table. The efficiencies (e^s) in Table [1] and 
Table [2] are defined as A^/ = eiiV*j_i, where i=l,2,3, Ni^ is the number of events selected 
after the ith cut out of A^^o* generated events. The expected number of events Ni in Table 
[1], where i=l,2,3, is obtained by multiplying the combined efficiency e(= €162. -.ei) with a C 

The expected number of signal events S (considering final states with both and /i^) is 
given by A.ePaC, where P is the product BR (see Eq.5) sufficient to yield S/\fB = 5. In our 
simulations e has been computed by generating {ti —>■ e^b) events only and we have assumed 
that the efficiency is the same for electrons and muons. Out of the three P's estimated 
in Table [2] the minimum {Pmin) is obtained by the combination of Ci and G2 only. This 
conclusion holds for other choices of mj^ (within the kinematic reach of Tevatron Run II) 
as well. We have tried a variety of additional selection criteria including b - tagging not 
shown in Table [1]. For example, we have tried various lower Pt cuts on the hardest jet. 
But the corresponding Pmin turns out to be weaker. On the other hand a lower Pt cut on 
second hardest jet yields a P sensitive to the LSP mass and introduces model dependence. 
We therefore conclude that the combination of Ci and C2 is the optimal one. 
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Backgrounds 


a{ pb) 






62 


N2 


63 


A^3 


tt 


3.73 


0.0179 


601.0 


0.7791 


468.3 


0.9327 


436.8 


w+w- 


9.51 


0.00623 


533.2 


0.2055 


109.5 


0.7188 


78.7 


wz 


1.16 


0.00541 


56.5 


0.2683 


14.9 


0.6831 


10.4 


WH 


0.11 


0.0093 


9.2 


0.3211 


3.0 


0.6286 


1.9 


bb 


2.82x10^ 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


tb +ib 


0.33 


0.00824 


24.47 


0.3398 


8.31 


0.8 


6.65 


Wbb 


11.6 


.00126 


131.5 


0.119 


15.7 


0.6667 


10.5 


Wuu 


73.5 


0.00103 


681.3 


0.2136 


145.5 


0.9545 


138.9 



Table 1: Leading order cross sections of the simulated backgrounds and the efficiencies of 
the cuts Ci, C2 and C3. 



Signal 


a{ pb) 




P 


£2 


P 




P 




0.41 


0.30473 


0.060 


0.9116 


0.037 


0.7702 


0.047 



Table 2: The signal cross section for mj^ = 180 GeV and efficiencies of the cuts Ci, C2 and 
C3 . We have computed P ( see Eq.5) after each cut from Eq.6 by requiring S/\fB = 5. 



In Table [2] P (Eq.5) is calculated by the formula: 



P = 



4>C(7(Fi?t)e, 



(6) 



where and denote the cross section and the combined efficiency (i.e., = 
N selected/ N generated — ^162) of the background of type b. Similarly e is the combined effi- 
ciency for the signal. The integrated luminosity C is taken to be 9 fb~^. 

We have also computed the backgrounds due to one valence quark and one sea quark. 
Some typical process and the corresponding number of background events (given in paren- 
theses) subject to the cuts Ci and C2 are: W"^dd^W-dd (12.4) and W^ud^W-ud (18.7). 
The corresponding crs have been calculated by applying nominal cuts of Pt > SGeV and 
\rj\ < 4.5 on the parton jets. It can be readily checked that these additional backgrounds 
hardly affect the estimated Pmm- 

We present in Fig.l Pmin as a function of m^^ for m~o= 120 GeV. In Fig. 2 the variation 
of Pmin with m~o is shown for mj^ —ISO GeV. It follows that Pmin is almost insensitive to 
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the LSP mass as claimed above. 




Figure 1: The model independent minimum value of the parameter P = BR(ii cXi) ^ 

BR(ti — i> Ifb), where /"*" = and iJ,^, observable at Tevatron Run II via 11 + jets + fjr 

channel as a function of mr . 

ti 

We present in Fig. 3 the distribution (unnormalized) of invariant mass of the highest Pt 
lepton and the highest Pt jet in the signal for mj^ —ISO GeV and m~o = 120 GeV. In spite 
of the presence of multiple jets in the signal, the combinatorial background does not obscure 
the peak at the input value of mj^ . In principle mj^ can be determined from this distribution 
and the lepton number violating nature of the underlying interaction can be established. 
Whether this peak will stand over the background depends on the actual value of P and mj^ . 
The issue of the observabihty of this peak at Tevatron, therefore, cannot be settled at the 
moment. 

In Fig. 4 we plot the distribution (unnormalized) of the invariant mass of the two lowest 
Pt jets in the signal, which comes most of the time from LSP decay. The mass of the LSP can 
be estimated in principle from the upper edge of this distribution. Again the combinatorial 
backgrounds obscure the edge a little bit but are not particularly severe. 

Throughout this section we assume BR(x5 — vhh)— 1.0. However, if this assumption 
is not strictly true, e.g., due to the presence of several rare decay modes of the LSP due to 
chargino-charge lepton or neutralino -neutrino mixing (see the next section) which we have 
neglected, then Pmin in Fig-1 should be interpreted as the minimum value of the product 
Px BR(x? ^ i^hl). 
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Figure 2: The variation of Pmin (see the caption of Fig.l) with m~o. 

3 Calculation of the product branching ratio in the 
model of m^y 

In this section we shall calculate P in some realistic models of neutrino mass constrained by 
the neutrino oscillation data and examine whether the predictions exceed Pmin estimated in 
the last section. Our main aim is to illustrate that the Run II data is sufficiently sensitive 
to probe these models and not to make an exhaustive study of all possible models. 

The collider signatures considered in the last section arise only in models with some non 
vanishing trilinear A^g^ type couplings. However, consistency with neutrino oscillation data 
require the introduction of more RPV parameters (bilinear superpotential terms, bilinear 
soft breaking terms etc). In fact the list of possible choices is quite long. It is expected that 
the constraints on the A' couplings in the most general model imposed by the u - oscillation 
data will be considerably weaker. On the other hand for large A' couplings the dilepton+dijet 
signal studied in earlier works [13] has a better discovery potential. Thus we have restricted 
ourselves to models with a minimal set of parameters capable of explaining the oscillation 
data with rather stringent constraints on the A' couplings. 

The benchmark scenario studied in this paper was proposed in ref. [U] where it was 
confronted with the then u - oscillation data. We work in a basis where the sneutrino vevs are 
zero. It is assumed that in this basis only three nonzero bilinear (/ij) and three trilinear (A -33) 
couplings, all defined at the weak scale, are numerically significant. In this framework the 
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Figure 3: The invariant mass distribution of the hardest jet and the hardest lepton in the 
signal for mj^ = 180 GeV. 

neutrino mass matrix receives contributions both at the tree and one loop level. It should be 
emphasized that the tree level mass matrix yields two massless neutrinos. Thus the interplay 
of the tree level and one loop mass matrices is essential for consistency with the oscillation 
data. 

The char gino- charge lepton, the neutralino - neutrino and other relevant mixing matrices 
in this basis may be found in [0]. In principle the diagonalization of these matrices may 
induce additional lepton number violating couplings which can affect the top squark and 
LSP decays considered in this paper. However, we shall show at the end of this section that 
the decays triggered by such induced couplings are highly suppressed either kinematically 
or dynamically in a wide variety of models. As a result the approximation that the decays 
of the top squark NLSP are driven by the A-33 couplings only is justified. The decay of the 
LSP requires a more careful handling and will be taken up at the end of this section. 

In [15] only the upper bounds on A-33 couplings as obtained from the neutrino oscillation 
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Figure 4: The invariant mass distribution of the lowest two Pt jets in the signal for m~o ~ 
120 GeV 

data were reported. It was shown in refs. [12] ([IS]) that the BR limits of the RPV decay of 
the ti - NLSP sensitive to Tevatron Run I data ( the MOBR that can be probed by Run II 
data) correspond to A's which are close to the above upper bounds. In ref. [17] the six RPV 
parameters in these models were randomly generated and the neutrino masses and mixing 
angles were computed for some well motivated choices of the RPC parameters. Compar- 
ing these with more recent oscillation data [28ja remarkably small allowed parameter space 
(APS) was obtained. It was also shown in [17] that there are six generic RPV scenarios 
consistent with the oscillation data ^28j. They are : 

a) fii < /i2,/i3: 

('^l) -^333 > -^133 ^ -^233 
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Figure 5: Number of points in the space allowed by the v oscillation data vs the parameter 
P(%) (Eq.5) in the gaugino like model with mq^ = 180 GeV^. For the choice of parameters 
and other details see text. 

('^2) A233 > A'i33 > A333 

b) fi2 < /ii,/i3: 

(^1) A'i33 ~ A233 >> A333 
(^2) A233 > A333 >> A'i33 

C) /i3 < fll,fi2- 

(Cl) A333 > A'i33 >> A233 

(C2) A333 > A233 >> A'i33 

Each scenario has its characteristic hierarchy among the three leptonic BRs of ti (see 
Eq.l). We focus on scenarios (61) and (62) which corresponds to relatively large BR for the 
decay channel in Eq.l with Zj = e or /x. However, in addition to the tree level and A' — A' 
loop contributions considered in |T5] we have included the contribution of the — A' loops 
to the z/ - mass matrix. An approximate form of the later can be found in|29[ [9] (see, e.g., 
Eq.35 of ref.[9]). The inclusion of the new contribution does not drastically alter the nature 
of RPV parameter space allowed by the z/-oscillation data. 

In addition to the RPV parameters the neutrino masses and mixing angles depends on 
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Figure 6: Number of points in the space allowed by the u oscillation data vs the parameter 
P(%) (Eq.5) in the mixed model with mj^ = 130 GeV. For the choice of parameters and 
other details see text. 

RFC parameters. In this paper we shall use the following popular assumptions to reduce 
the number of free parameters in the RFC sector: i) At the weak scale the soft breaking 
mass squared parameters of the L and R-type squarks belonging to the third generation are 
assumed to be the same( the other squark masses are not relevant for computing neutrino 
masses and mixing angles in this model), ii) We shall also use the relation M2 ~ 2 Mi at 
the weak scale as is the case in models with a unified gaugino mass at Mq- Here Mi and 
M2 are respectively the soft breaking masses of the U(l) and SU(2) gauginos respectively. 

The tree level neutrino mass matrix and,hence, the predicted neutrino masses depends 
on the parameters of the gaugino sector (through the parameter C fTEl [T7]). They are M2, 
Ml , /i (the higgsino mass parameter) and tan /3 = V2/V1, where vi and V2 are the vacuum 
expectation values (vevs) for the down type and the up type neutral higgs bosons respectively. 
We remind the reader that for relatively large tan j3s the loop decay overwhelms the RFV 
decay. We have, therefore, restricted ourselves to tanP = 5-8. 

It is also convenient to classify various models of the RFC sector according to the relative 
magnitude of M2 and /i. If Mi < M2 <^ /i, then the lighter chargino ( xf ), the LSF (xi) ci-nd 
the second lightest neutralino (X2) dominantly gauginos. Such models are referred to as 
the gaugino-like model. On the other hand in the mixed model( Mi< M2~ /i), xf and X2 
are admixtures of gauginos and higgsinos. In both the cases, however, Xi is almost a bino. 
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There are models with Mi,M2 ^ /i in which xf^ Xi ^^'^ X2 higgsino - hke and all have 
approximately the same mass ( ~ yu). It is difficult to accommodate the top squark NLSP 
in such models without fine adjustments of the parameters. Thus the LSP decay seems to 
be the only viable collider signature. 

One can also construct models wino or higgsino dominated LSPs. However, the ti-NLSP 
scenario cannot be naturally accommodated in these frameworks for reasons similar to the 
one in the last paragraph. 

The one loop mass matrix, on the other hand, depends on the sbottom sector (through 
the parameter K2 p!5|[T7]). This parameter decreases for higher values of the common squark 
mass for the third generation. From the structure of the mass matrix it then appears that 
for fixed C, identical neutrino masses and mixing angles can be obtained for higher values 
of the trilinear couplings if K2 is decreased. Thus at the first sight it seems that arbitrarily 
large width of the RPV decays may be accommodated for any given neutrino data. This, 
however, is not correct because of the complicated dependence of the RPV and loop decay 
BRs of ti on the RPC parameters and certain theoretical constraints. The common squark 
mass cannot be increased arbitrarily without violating the top squark NLSP condition. Of 
course larger values of the trilinear soft breaking term At may restore the NLSP condition. 
But larger values of At tend to develop a charge colour breaking( CCB ) minimum of the 
scalar potential ^30j. Finally the pseudo scalar higgs mass parameter Ma can be increased 
to satisfy the CCB condition. But as noted earlier that would enhance the loop decay 
width as well and suppress the BRs of the RPV decay modes. 

We first examine scenario(62) with the following set of RPC parameters: 
A) Ml = 100.0, M2 = 200.0, = 320.0, tan /? = 7.0, A = 938.0, ^ = 300.0, M^- = 400.0, M,- 
(common slepton mass ) = 350.0 and Ma = 600.0, 

where all masses and mass parameters are in GeV |. We note in passing that common 
slepton mass does not enter into the calculation of the BRs and any choice which preserves 
of ti - NLSP condition serves the purpose. For the parameters chosen the loop(4-body) 
decay BR is ^ 50% (34%). Yet the bulk of the APS gives P > P™„. 

In this gaugino like model mj^ = 180 GeV and m~o=100 GeV. The corresponding model 
independent Pmin is 0.037(see Fig.l). 

We then randomly generate 10^ sets of the six RPV parameters under consideration and 
count the sets allowed by the oscillation data. For each point of the APS we compute BRs 
of the three modes in Eqs.1,2 and 3 and get the corresponding P. 

In Fig. 5 we present a histogram of the number of sets allowed by the oscillation data 
vs P in the above gaugino like model. It is found that for most of the APS P > Pmin- It 
may be recalled that in ref. [TT] it was pointed out that in this model the loop decay BR is 
much larger than the total RPV decay BR of ti even for modest values of tan (3. Thus the 
inclusion of the fi — X' loop does not change our conclusion drastically. 

In Fig. 6 we present a similar histogram in a mixed model. Here the LSP has a relatively 

^For computing the loop decay BR the mass of the charm squark is required. We assume it to be equal 
to the third generation common squark mass. 
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large higgsino component. Thus the additional decay modes discussed in the introduction 
(see after Eq.4) may open up and reduce the combined BR of the modes in Eq. 4. Hence 
we have restricted ourselves to m-o < mw- 

We have considered the {hi) scenario for the RPV parameters. The RPC parameters 
corresponding to this mixed model are chosen to be: 

B) Ml = 78.0, M2 = 170.0, /i = 180.0, tan /5 = 8.0, A = 890.0, A;, = 1000.0, = 
375.0, Ml = 350.0, Ma = 300.0, 

where all mass and mass parameters are in GeV. With the above choice of parameters mj_^ = 
130 GeV and m~o. Since, as discussed in section 2, Pmin is highly insensitive to m~o one can 
still use the estimates of figure 2. In this case the loop(4-body) decay BR is ~ 85%(6.5%). 
Yet the entire APS gives P > Pmin- 

We next compare and contrast the signals considered in this paper and the one discussed 
in [13]. In the latter work assuming both the tis produced at Tevatron Run II decay via 
the tl e^d channel the MOBR of this channel was estimated to be 20%. The Drell-Yan 
process turned out to be the dominant source of background. It was further noted that if 
b -tagging is employed the signal will be essentially background free. Using the efficiencies 
as given in Table [2] of ref.[T3] and including a b-tagging efficiency of 50% the MOBR for 
the channel ti e^h corresponding to ten signal events was estimated to be roughly the 
same. We define the parameter BR(ti — »• e^b) + BR( ti /i^&)= BR(e + /x). Following 
the procedure of ref. briefly sketched above, we then estimate the minimum observable 
value of this parameter (= MOBR(e + yu)) for £ = 9 fb~^. The results are given in Table [3] 





a{ pb) 


e 


MOBR(e + ^) 


100 


13.1 


0.0194 


0.076 


140 


2.1 


0.0933 


0.087 


180 


0.41 


0.2278 


0.126 


220 


0.12 


0.3073 


0.20 



Table 3: The minimum observable value of BR(ti e~^b) +(ti — fi'^b) at Tevatron Run II 
via the dilepton dijet channel for different mj^. 

We flnd that in any gaugino like model with tan (5 > 6 the computed BR(e + /i)s turn out 
to be smaller than the MOBR in Table [3] practically over the entire APS . For example, in 
the gaugino like model considered above with tan /3 = 7, BR(e) ~ and BR(/i) varies from 
0.06 to 0.11 (except for a few solutions which corresponds to BR(/i) around 0.025). Thus 
BR(e + /i) is indeed much smaller than the MOBR in Table [3] for mj^ = 180 GeV. However, 
most of the points in the APS yield P larger than Pmm(Fig.5). 

For the mixed model considered above the entire APS yields P > Pmin- On the contrary 
BR(e + /i) is still below the MOBR since the loop and 4-body decay BRs are ~ 85% and 
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~ 6.5% respectively. These examples illustrate that the signal in [13] and the one considered 
in this paper are indeed complimentary. 

As noted at the beginning of this section some RPC interactions can induce new lepton 
number violating interactions of ti and Xi due to chargino- lepton or neutralino-neutrino 
mixing. These induced interactions can in principle affect the decays of the ti-NLSP or the 
LSP considered in section 2. For example the RPC vertex ti - h - ( or H^) may lead 
to additional lepton number violating couplings of the ti due to chargino-lepton mixings. 
Similar induced couplings may arise from ti - t - Wsl or B or H^) coupling due to mixing in 
the neutralino-neutrino sector ( these couplings may be relevant only if the ti is significantly 
heavier than the t). The mixing factors which would suppress the induced couphngs can be 
estimated from the 5x5(7x7) chargino - lepton ( neutralino-neutrino) mass matrix. The 
estimated value is O{fii/100GeV) where 100 GeV is the typical magnitude of an element 
of the 2x2 chargino block or the 4x4 neutralino block of the above matrices. Since the 
largest fii allowed by the oscillation data is C(10~^) GeV the mixing factors are estimated 
to be O{10~^) or smaller. Moreover, the induced couplings will be additionally suppressed 
by gauge or Yukawa couplings. On the other hand the smallest A^gg coupling contributing 
to ti decay consistent with oscillation data is O{10~^). Thus the ti-NLSP BRs computed 
by considering \[^^ driven decays only are quite reliable. The rough estimates presented 
here would be substantiated below by results obtained by numerically diagonalizing the 
chargino-lepton and neutralino-neutrino mass matrices. 

The LSP decays require more careful analysis. The main decay mode considered in this 
paper (Eq 4) is a three body decay. On the other hand lepton number violating two body 
decays of the LSP (the decays Xi ~^ and Xi ~^ /"""W^ are examples) can be induced by 
the Xi - - W~ or Xi - ^3 - Z vertices. However, when the LSP is almost a pure bino, 
which is the case in the gaugino model, the original RPC couplings are highly suppressed. 
In addition suppression by the mixing factors discussed in the last paragraph will come into 
play. Consequently the BRs of the lepton number violating 2 body decays of the LSP are 
0(10.0%) (see below for numerical results). Thus the decay in Eq. 4 indeed occur with 
almost 100 % BR in this model. 

In more general models with the LSP having significant higgsino components (e.g, in the 
mixed model) our assumption regarding the LSP decay is valid if m~o < m\y. For heavier 
LSPs the parameter Pmin in section. 2 should be interpreted as the minimum value of the 
product P xBR{xi ^ibb) observable at Run II. 

We now present some numerical results. We numerically diagonalize the mass matrices 
in the chargino-lepton or neutralino-neutrino sector for all combinations of RPV parameters 
allowed by the oscillation data. For the parameter set A) (the gaugino like model) we find 
that the maximum amplitude for finding a charge lepton mass eigenstate in a or is 
3.3 xlO"^. The corresponding A233, responsible for the ii or LSP decay signal, is 8^ xlO^. 
Similarly the maximum amplitude for finding a neutrino mass eigenstate in a B, or 
is 3.6 xlO~^. The corresponding A233 is 8.6 xlO~^. Using the induced couplings in this 
scenario as given above, the mixings in the RPC chargino and neutralino mass matrices and 
the widths of the modes {Xi ~^ ^x'j) and {xa ~^ W'^lj'^) given, e.g., in [21] Eq.l5 to Eq.21, 
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we find that BR(x° Zuj) = 13.0% and BR(x- ^ W+lj) = 1.0%. 



4 Conclusion 

In conclusion we reiterate that the 11+ jets(> 3)+^ signal arising from RPC and RPV 
decays of ti - tl pairs produced at the Tevatron is a promising channel for probing a class 
of RPV models of neutrino mass. For a set of kinematical cuts suitably optimized in this 
paper the size of this signal is essentially controlled by the production cross section of the ti 
- tl pair as given by QCD and the parameter P (see Eq.5). Using Monte Carlo simulations 
we have obtained model independent estimates of Pmin for different mj^s corresponding to 
an observable signal for an integrated luminosity of 9 fb~^ (see Fig.l). The efficiencies of 
the cuts are entirely controlled by the kinematics of the lepton and the jet with highest 
which most of the time come directly from the RPV decay of the ti. The size of the signal 
and, consequently, the estimated Pmin are practically independent of m~o (Fig.2). 

We have also noted that in spite of the combinatorial backgrounds, the invariant mass 
distribution of the above lepton-jet pair shows a peak at m^^(see Fig. 3). This peak, if discov- 
ered, will clearly establish the lepton number violating nature of the underlying interaction. 
Similarly the end point of the invariants mass distribution of the two lowest jets (Fig. 4) 
which principally arise from LSP decays ( Eq.4) may determine the LSP mass. Whether 
these salient features of the distributions will be obscured by the full SM background in the 
Tevatron data depends on the actual values of P and mj^ . 

This signal may turn out to be the main discovery channel even if the loop decay (Eq.2) 
of the ti-NLSP (followed by the LSP decay) strongly dominates over its RPV decay (Eq.l). 
On the other hand if the RPV decay mode overwhelms the loop decay then the dilepton 
+ dijet channel studied in [HI [121 [E] may provide a better signal. Finally if the data 
establish a competition between the two modes that would also be highly indicative of an 
underlying RPV model of neutrino mass. It may be recalled that in these models the neutrino 
oscillation data requires the A^gg couplings to be highly suppressed. As a result the 2-body 
RPV decays have widths comparable to the competing RPC decays which occur in higher 
orders of perturbation theory if ti is the NLSP. 

The prospect of discovering the RPV model considered in this paper will be better if 
the RPV decays of ti into final states with r+b can also be probed at the Tevatron. In fact 
scrutinizing these models in the light of the oscillation data reveal that ti T~^h is indeed 
the most dominant decay mode over a large region of the APS |I7]. A model independent 
estimate of Pmin for this channel is, therefore, very important for a complete probe of this 
model. 

Our computation of P in specific benchmark models [151 [IB] of neutrino oscillations 
establish that the APS filtered out by the oscillation data contain many points with P> Pmin- 
This happens both in the gaugino like( Fig.5) and in the mixed model (Fig. 6). There are 
also regions in the APS where the dilepton + dijet signal proposed in [13] is unobservable 
but the signal proposed in this paper stands over the background. Thus the two signals are 
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indeed complementary. 

Since in the RPV MSSM leptons are baryons must be treated differently there is a basic 
incompatibility between this model and a typical GUT in which quarks and leptons are placed 
in the same multiplet. Thus RPV terms in the superpotential tends to violate both baryon 
and lepton number conservation and lead to catastrophic proton decays. Thus the task of a 
model builder is to remove such terms by introducing appropriate discrete symmetries [52] . 
Yet the RPV MSSM can be accommodated in the framework of a GUT. RPV interactions 
may, for example, be induced at the GUT scale by higher dimensional non-renormalizable 
operators which reduce to either baryon number or lepton number violating interactions 
when the GUT breaks down to the SM and certain heavy scalar fields develop vevs [33] • 

In GUT models neutrino masses can be generated at the weak scale in a variety of 
interesting ways. The set of input RPV parameters at Mq need not be identical to the set 
appearing in the neutrino mass matrix at the weak scale. In fact the former set may have 
smaller number of parameters than the latter set. For example, one may start at Mq with 
three relatively large trilinear couplings different from the A^ggS required by the neutrino 
sector. Renormalization group evolution [51] and flavour violation inevitably present in any 
model due to the CKM mixing would then induce the A^gg and the /Xj parameters at the 
weak scale [35]. However, other RPV parameters may also be generated leading to a more 
complicated v- mass matrix. The relatively large input couplings may then lead to a rich 
low energy phenomenology [35] in addition to ti and LSP decays. In this paper, however, 
we have not considered the origin of the weak scale parameters and have restricted ourselves 
to the signatures of the A^gg couplings only. 

The signal discussed in this paper and the one in ref. [T3] will certainly have much larger 
sizes at the LHC. But at higher energies many other sparticles may be produced as well. 
Thus one has to isolate the signal not only from the SM background but also from the SUSY 
background. On the other hand since the ti - NLSP may very well be the only strongly 
interacting sparticle within the kinematic reach of the Tevatron. As a result these signals 
may be observed in a relatively clean environment. 
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